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ABSTRACT 



Context. The lack of knowledge of photospheric parameters and the level of chromospheric activity in young low-mass pre-main 
sequence stars introduces uncertainties when measuring mass accretion rates in accreting (Class II) Young Stellar Objects. A detailed 
investigation of the effect of chromospheric emission on the estimates of mass accretion rate in young low-mass stars is still missing. 
This can be undertaken using samples of young diskless (Class III) K and M-type stars. 

Aims. Our goal is to measure the chromospheric activity of Class III pre main sequence stars to determine its effect on the estimates 
of accretion luminosity (L acc ) and mass accretion rate (M acc ) in young stellar objects with disks. 

Methods. Using VLT/X-Shooter spectra we have analyzed a sample of 24 non-accreting young stellar objects of spectral type between 
K5 and M9.5. We identify the main emission lines normally used as tracers of accretion in Class II objects, and we determine their 
fluxes in order to estimate the contribution of the chromospheric activity to the line luminosity. 

Results. We have used the relationships between line luminosity and accretion luminosity derived in the literature for Class II objects 
to evaluate the impact of chromospheric activity on the accretion rate measurements. We find that the typical chromospheric activity 
would bias the derived accretion luminosity by L acclloisi ,< 10~ 3 L o , with a strong dependence with the T cff of the objects. The noise on 

M acc depends on stellar mass and age, and the typical values of log(M acc-no i sc ) range between 9. 2 to -11 .6M Q /yr. 

Conclusions. Values of L acc < 10~ 3 L o obtained in accreting low-mass pre main sequence stars through line luminosity should be 
treated with caution as the line emission may be dominated by the contribution of chromospheric activity. 

Key words. Stars: pre main sequence — Stars: low-mass - Stars: activity - accretion 
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1. Introduction 

Circumstellar disks are formed as a natural consequence of an- 
gular momentum c onservation du ring gravitational collapse of 
cloud cores (e.g. I Shu et all Il987l) . In the early phases of star 
formation, the disk allows for the dissipation of angular mo- 
mentum channeling the accretion of material from the infalling 
envelope onto the central young stellar object (YSO). At later 
stages, when the envelope is dissipated, planetary systems form 
in the disk, while the star-disk interaction continues through the 
inner disk and the stellar magnetosphere. This phenomenon con- 
strains the final stellar mass build-up (e.g. Hartm ann et al.L fl998) 
and its typical timescales are connected with the timescales 
on which disks dissipate and plan e t formation occurs (e.g . 
Hernandez etall l2007t iFedele et all 12011 IWilliams & Ciezal 
201 lh . 



Accretion can be observed using typical signatures in the 
spectra of YSOs, such as t he continuum excess in the blue part of 
the visible spectrum (e.g. iGullbring et all Il998[) and the promt 
nent op tical and in fra red em i ssion lines (e~g. IMuzerolle et al 
1998 : INatta et all l2004t iHerczeg & H illenbrand. 2008; 
iRigliaco et all 120121) . The measurements used to determine 
the bolometric accretion luminosity (L acc ) are either direct 
or indirect. Direct measurements are obtained by measuring 
the emission in excess of the photospheric one in the Balmer 
and Paschen continua and adopting a model to correct for the 
emission at the wavelengths not covered by the observations, 
and which o r iginates mostly below the U - band threshold (e.; 
Valenti et al 



1998 



* Based on observations collected in the programs 084.C-0269 
085.C-0238, 086.C-0173, 087.C-0244, 089.C-0143 at the European |2()()5URigliaco et alll2()Ti1h 
Organisation for Astronomical Research in the Southern Hemisphere 
(Chile). 



1991 IGullbring et all [19981 ICalvet & Gul lbrins 
Herczeg & Hillenbrand 120081; [Rigliaco et all I20K 
or by line profile modelling (e.g. IMuz erolle et al. . 1998) 
Indirect measurements are obtained using empirical correla- 
tions between emission line luminosity (Ln ne ) and L acc (e.g . 



Muzerolle et all 119981; INatta et all 12004 12006; Mohantv etall 
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Measurements of mass accretion rate (M acc ) are subjected to 
many uncertainties, because this quantity depends on L acc and 
on the mass-to-radius (M*/R*) ratio. This ratio is normally de- 
termined from the position of the object on the HR diagram and a 
set of evolutionary models. Uncertainties on spectral type (SpT) 
of the objects affect the determination of effective temperature 
(T e ff) while those on the extinction (Ay) and the distance mainly 
affect the estimate of stellar luminosity (L»). It is not trivial to 
determine those parameters in accreting stars because of accre- 
tion shocks on the s tellar surface produc i ng ve iling in the pho- 
tospheric lines (e.g. ICalvet & Gullbringl fl998) and modifying 
the photometric colors (e.g. lDa Rio et al.L 120101) . Moreover, the 
derivation of Li; ne , from which L acc is determined, is affected 
by another stell ar property, namely the chromospheric activ - 
ity of the YSOs dHoudebine et all \l99^L iFranchini et all Il998l) . 
Chromospheric line emission is usually small compared to the 
accretion-powered emission, but it can become i mportant when 
accre tion decreases at later evolutionary stages (lngleby et al., 
1201 B and in lower m ass stars, where accretion rates are smaller 
dRigliaco et all 120121) . Therefore, this is an important source of 
uncertainty that has not been investigated in detail so far. 

Part of the INAF consortium's Guaranteed Time 
Observations (GTO) of X-Shooter, a broad-band, medium- 
resolution, high-sensitivity spectrograph mounted on the 
ESO/VLT, has been allocated for star formation studies, in 
particular to investigate accretion, outflows, and chromospheric 
emissio n in low-mass Clas s II young stellar and sub-stellar 
objects dAlcala et all 1201 ll) . The targets observed during the 
GTO were chosen in nearby (d < 500 pc) star forming regions 
with low extinction, and with many very low-mass (VLM) YSO 
(M» < 0.2 M G ). Generally, those YSOs for which measurements 
in many photometric bands were available, both in the IR (e.g. 
Hernande z et all 120071; iMerin et all 120081) an d in the visible 
part o f the spectrum (e.g. [Merinetal., 2008; Rigl iaco et all 
1201 lh were selected. In order to derive L acc of a given Class 
II YSO a Class III template of the same SpT as the Class II is 
needed. Therefore, during this GTO survey 24 Class III targets 
in the range K5-M9.5 have been observed, providing the first 
broad-band grid of template spectra for low-mass stars and 
brown dwarfs (BDs). Since these spectra have a very large 
wavelength range (~350 - 2500 nm) covering part of the UV 
spectrum (UVB), the whole Visible (VIS), and the Near Infrared 
(NIR), this sample allows us to: determine the stellar parameters 
of the targets, derive the chromospheric emission line fluxes 
and luminosities and, hence, determine the implications of 
chromospheric emission on the indirect accretion estimates in 
Class II objects. 

The paper is structured as follows. In Sect. |2]we discuss the 
sample selection, the observation strategy and the data reduc- 
tion procedure. In Sect. [3] we describe how SpTs of our targets 
have been determined, while in Sect. [4] we derive their main stel- 
lar parameters. In Sect. [5] we identify the main lines present in 
the spectra and derive their intensities. In Sect.[6]we discuss the 
implications of the line luminosity found for studies of M acc in 
Class II YSOs. Finally, in Sect.Qwe summarize our conclusions. 



2. Sample, observations and data reduction 

2.1. Sample 

Among the objects observed in the GTO survey, we selected 
only those that have been classified as Class III objects using 
Spitzer photometric data. The sample comprises Class III YSOs 
in the cr Orionis, Lupus III and TW Hya associations. In the 



end, the number of targets is 24; 13 objects are members of the 
TW Hya association, 6 of the Lupus III cloud and 5 of the <x 
Orionis region. Their SpTs range between K5 and M9.5 (see 
Sect. E). Three BDs, namely Par-Lup3-1, TWA26, and TWA29, 
are included in our sample, with SpT M6.5, M9, and M9.5, re- 
spectively. Data available from the literature for these objects are 
reported in Table Q] 

Six YSOs in our sample are components of three known 
wide visual binary systems. In all cases we were able to resolve 
them, given that their separations are always larger than 6". 

2.2. Observations 

All the observations have been made in the slit nodding mode, in 
order to achieve a good sky subtraction. Different exposure times 
and slit dimensions were used for different targets, in order to 
have enough S/N and to avoid saturation. The readout mode used 
in all the observations was "100,1x1, hg", while the resolution 
of our spectra is R = 9100, 5100, and 3300 in the UVB arm 
for slit 0.5", 1.0", and 1.6", respectively; R = 17400, 8800, and 
5400 in the VIS arm for slit 0.4", 0.9", and 1.5", respectively; 
R = 11300, 5600, and 3500 in the NIR arm for slit 0.4", 0.9", 
and 1.5", respectively. We report in Table [2] the details of all 
observations performed for this work. 

2.3. Data reduction 

The data reduction has been done u s ing tw o versions of the 
X-Shooter pipeline (Modigl iani et all 1201 Oh . run through the 
EsoRex tool, according to the period in which the data were ac- 
quired: version 1.0.0 has been used for the data of December 
2009 and May 2010, while for data gathered in January 2011, 
April 201 1, and April 2012, version 1.3.7 has been used. The two 
versions lead to results that are very similar. The reduction has 
been done independently for each spectrograph arm. This takes 
into account also the flexure compensation and the instrumental 
profile. We have used the pipeline recipe XSh_SCired_slit_nod 
which includes bias and flat-field correction, wavelength calibra- 
tion, order tracing and merging and flux calibration. Regarding 
the last point, by comparison of the response function of differ- 
ent flux standards observed during the same night, we estimate 
an intrinsic error on the flux calibration of less than 5%. Given 
that some observations have been done with poor weather condi- 
tions (seeing ~3.5") or with narrow-slits, we have then checked 
the flux calibration of each object using the available photomet- 
ric data, usually in the U, B, V, R, I, J, H, K bands, as reported in 
Table Q] We have verified that all the spectra match well the 
photometric spectral energy distribution (SED) and adjusted the 
flux-calibrated spectra to match the photometric flux. Binaries 
have been reduced in stare mode. Telluric removal has been done 
using standard telluric spectra obtained in similar conditions of 
airmass and instrumental set-up of the target observations. This 
correction has been accomplished with the IRAFQ task telluric, 
using spectra of telluric standards from which photospheric lines 
were removed using a multigaussian fitting. The correction is 
very good at all wavelengths, with only two regions in the NIR 
arm (AA 1330-1550 nm, AA 1780-2080 nm) where the telluric 
absorption bands sa turate. More detail about the reduction will 
be reported in lAlcala et al.l (12013. in prep.l) . 



1 IRAF is distributed by National Optical Astronomy Observatories, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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Table 1. Known parameters from the literature. 
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Table 2. Details of the observations. 
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16-17 May 2010 


Szl22 


1.0" 


0.9" 


0.9" 


600s 


18 Apr 2012 


TWA9B 


0.5" 


0.4" 


0.4" 


800s 


16-17 May 2010 


TWA15B 


0.5" 


0.4" 


0.4" 


600s 


16-17 May 2010 


TWA7 


0.5" 


0.4" 


0.4" 


100s 


12 Jan 2011 


TWA15A 


0.5" 


0.4" 


0.4" 


600s 


16-17 May 2010 


Szl21 


1.0" 


0.9" 


0.9" 


500s 


18 Apr 2012 


Sz94 


1.0" 


0.9" 


0.9" 


600s 


16-17 May 2010 


S0797 


1.0" 


0.9" 


0.9" 


2400s 


23-24 Dec 2009 


SU641 


1.0 


0.9 


n o" 
0.9 


360US 


23-24 Dec 21)09 


Par-Lup3-2 


1.0" 


0.9" 


0.9" 


1200s 


16-17 May 2010 


S0925 


1.0" 


0.9" 


0.9" 


3600s 


21-22 Dec 2009 


S0999 


1.0" 


0.9" 


0.9" 


2400s 


24-25 Dec 2009 


Szl07 


1.0" 


0.9" 


0.9" 


600s 


22 Apr 2011 


Par-Lup3-1 


1.0" 


0.9" 


0.9" 


600s 


16-17 May 2010 


TWA26 


1.0" 


0.9" 


0.9" 


3600s 


22 Mar 2010 


TWA29 


1.6" 


1.5" 


1.5" 


3600s 


22 Mar 2010 



The fully reduced, flux- and wavelength-calibrated spectra 
are available on the ESO Archival 

3. Spectral type classification 

A careful SpT classification of the sample is important in order 
to provide correct templates for accretion estimates of Class II 
YSOs. Moreover, the procedure used to derive the SpT of Class 
II and Class III YSOs should be as homogeneous as possible. 
In this section we describe two different methods to derive SpT 
for these objects. Firstly, we use the depth of various molecular 
bands in the VIS part of the spectrum. Then, we describe the 
second method, which consists of using spectral indices in the 
VIS and in the NIR part of the spectrum. These provide us a 
reliable, fast and reddening free method to determine SpT for 
large samples of YSOs. 

3. 1 . Spectral typing from depth of molecular bands 

For the SpT classification of the objects, we use the analysis 
of the depth of several molecu lar bands in the spectral region 
between 580 nm and 900 nm (iLuhmanl I20041 lAllen & Stroml 
ll995_lHenrvet al.lll994l) . This region includes various TiO ( AA 
584.7-605.8, 608-639, 655.1-685.2, 705.3-727, 765-785, 820.6- 
856.9, 885.9-895 nm), VO (AA 735-755, 785-795, 850-865 nm) 
and CaH (AA 675-705 nm) absorption bands, and a few photo- 
spheric lines (the Call IR triplet at AA 849.8, 854.2, 866.2 nm, 
the Nal doublet at A 589.0 and 589.6 nm, the Cal at A 616.2 nm, 
a blend of several lines of Ball, Fel and Cal at A 649.7 nm, the 
Mgl at A 880.7 nm, and the Nal and KI doublets at AA 818.3 nm 
and 819.5 and AA 766.5 nm and 769.9, respectively. 

In Fig.[TJ |2]and [3]we show the VIS spectra of the objects in 
the wavelength range between 580 and 900 nm. All the spectra 
are normalized at 750 nm and, for the sake of clarity, smoothed 

2 http://www.eso.org/sci/observing/phase3/data_releases.html 




K3 K5 K7 Ml M3 M5 M7 M9 
SpT 



Fig. 4. Distribution of spectral types of the Class III YSOs dis- 
cussed in this work. Each bin corresponds to one spectral sub- 
class. 



to a resolution of R~2500 at 750 nm and vertically shifted. The 
depth of the molecular features increases with SpT almost mono- 
tonically, and, comparing the spectra of the targets using together 
different wavelength subranges and different molecular bands, 
we can robustly assign a SpT to our objects thanks to the differ- 
ences in the depth of the bands, with uncertainties estimated to 
be 0.5 subclasses. 

With this spectral typing procedure, we classify all the ob- 
jects in our sample with SpT earlier than M8. The agreement 
between the SpTs derived here and those in the literature is 
good, with only three cases where the difference is 2 spec- 
tral subclasses. We assume the SpT availabl e from the litera- 
ture for the 2 YSOs wit h SpT later than M8 dReid et all [2008; 
Kirk patrick et all 120081) . because the classification by compari- 
son of molecular bands depth with other spectra in the sample is 
not possible due to the fact that our sample has a gap between 
M6.5 and M9. We can only confirm that these objects have a SpT 
later than the other targets in the sample and that their SpT dif- 
fer by 0.5 subclass. The SpT obtained here are listed in Table [3] 
and in the first column of Table |4] and they are used for the rest 
of the analysis. The distribution of SpT for the sample is shown 
in Fig. [4] The range of the M-type is almost entirely covered, 
providing a good sample for the goals of this paper and a solid 
library of templates that can be used for L acc estimates of Class 
II YSOs. 

In Appendix [C] the NIR and UVB spectra of all the sources 
are shown. Also in these cases a trend with SpT can be seen. 

3.2. Spectral indices for M3-M8 stars 

Spectral indices provide a fast method to determine SpT for large 
samples of obj ects. Here we test some of these indices for M- 
type objects. Riddi ck et alJ ([2007) tested and calibrated various 
spectral indices in the VIS part of the spectrum for pre main se- 
quence (PMS) stars with SpT from M0.5 to M9. In fact, they 
suggest the use of some reliable spectral indices valid in the 
range M3-M8. We find that the best SpT classification can be 
achieved by combining together results obtained using the set 
of indices that we report in Table [5] We proceed as follows. For 



4 



Manara et al.: Impact of chromospheric emission on accretion rate estimate of young stellar objects 




Wavelength (nm) 



Fig. 1. Spectra of Class III YSOs with SpT earlier than M3 in the wavelength region where the spectral classification has been 
carried out (see text for details). All the spectra are normalized at 750 nm and offset in the vertical direction by 0.5 for clarity. The 
spectra are also smoothed to the resolution of 2500 at 750 nm to make easier the identification of the molecular features. 

each object we calculate the SpT with all these indices, and then the final results obtained with these indices in the second col- 
assign the mean SpT using those results that are in the nominal umn (VIS Jnd) of Table [4] Comparing these results with those 
range of validity of each index. Typical dispersions of the SpT derived in the previous section, we report an agreement within 
derived with each index are less than half a subclass. We report one subclass for all the objects in the range M3-M8, as expected. 
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Fig. 2. Same as Fig.Q] but for SpTs between M3 and M5. 



However, there are 4 YSOs classified M1-M2 from the depth of 
molecular bands that would be classified M3 from the values of 
the Riddick's indices. This suggests that spectral classification 
M3 obtained with spectral indices could, in fact, be earlier by 
more than one subclass. 



There are also indices based on features in the NIR par t of the 
spectrum (e.g Testi et al., 2001; Testi, 2009; Allers et al 1 120071 
iRoias-Avala et al.Ll2012l) . In AppendixlBlwe compare the results 
obtained with these NIR indices with the SpT derived in the pre- 
vious section, in order to confirm the validity of some of these 
indices for the classification of M-type YSOs. 
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Wavelength (nm) 



Fig. 3. Same as Fig. Q] but for SpTs later than M5. 



4. Stellar parameters 

To estimate the mass (M»), radius (R*) and age for each target 
we compare the photosp heric parameters (L, ,T P ,ff) with the theo- 
retical predictions of the Bar affe et al.l d!998l) PMS evolutionary 
tracks. We derive the T e ff of each star from its SpT using the 



iLuhman et al.1 (120031) SpT-T e ff scale, while the procedure to esti- 
mate L„ is described in the next paragraph. We then use (L,,T e ff) 
to place the stars on the HR diagram and estimate (M„, age) in- 
terpolating the theoretical evolutionary tracks. These parameters 
are reported in Table [3] Finally, we obtain R* from L, and T e ff. 
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Table 3. Stellar parameters derived for the objects in our sample 



Name 


SpT 


Teff [K] 


log(L,/L ) 


M.[M ] 


<log(L„„..„ , Jt ,/L Q )> 


TWA9A 


K5 


4350 


-0.61 


0.81 


-3.38 


S0879 


K7 


4060 


-0.29 


1.07 


-2.88 


TWA6 


K7 


4060 


-0.96 


0.66 


-3.37 


TWA25 


MO 


3850 


-0.61 


0.84 


-3.04 


TWA 14 


M0.5 


3780 


-0.83 


0.73 


-3.23 


TWA13B 


Ml 


3705 


-0.70 


0.68 


-3.32 


TWA13A 


Ml 


3705 


-0.61 


0.70 


-2.85 


TWA2A 


M2 


3560 


-0.48 


0.55 


-3.36 


Szl22 


M2 


3560 


-0.60 


0.54 


-2.57 


TWA9B 


M3 


3415 


-1.17 


0.37 


-3.93 


TWA15B 


M3 


3415 


-0.96 


0.37 


-3.45 


TWA7 


M3 


3415 


-1.14 


0.37 


-3.84 


TWA15A 


M3.5 


3340 


-0.95 


0.30 


-3.23 


Szl21 


M4 


3270 


-0.34 


0.37 


-3.25 


Sz94 


M4 


3270 


-0.76 


0.28 


-3.38 


S0797 


M4.5 


3200 


-1.26 


0.19 


-4.27 


S0641 


M5 


3125 


-1.53 


0.12 


-4.51 


Par-Lup3-2 


M5 


3125 


-0.75 


0.18 


-3.89 


S0925 


M5.5 


3060 


-1.59 


0.10 


-4.65 


S0999 


M5.5 


3060 


-1.28 


0.13 


-4.30 


Szl07 


M5.5 


3060 


-0.79 


0.16 


-3.69 


Par-Lup3-1 


M6.5 


2935 


-1.18 


0.10 


-4.74 


TWA26 


M9 


2400 


-2.70 


0.02 


-6.54 


TWA29 


M9.5 


2330 


-2.81 


0.02 


-6.64 



Notes. The spectral type-T c ff relation is adopted from lLuhman et all (12003) for M-type objects and from Kenvon & Hartmann ( 1995) for K-type 
objects. 



4.1. Stellar luminosity 

Given the large wavelength coverage (~300-2500 nm) of the X- 
Shooter spectra, for objects with 2300 < T eff < 4400 K only 
a small percentage of the stellar flux (< 10-30%) arises from 
spectral regions outside the X-Shooter spectral range. We use, 
therefore, the following procedure to estimate the total flux of 
our objects: first, we integrate the whole X-Shooter spectrum 
from 350 nm to 2450 nm, excluding the last 50 nm of spectra 
on each side, that are very noisy, and the regions in the NIR 
between the J, H and A"-bands (AA 1330-1550 nm, AA 1780-2080 
nm, see Sect. [2]), where we linearly interpolate across the telluric 
abso rption region s . We t hen use the BT-Settl synthetic spectra 
from Alla rd et all Q201 1 ) with the same T e ff as our targets (see 
Table [3]), assuming logg=4.0, typical of low-mass YSOs, and 
normalized with our spectra at 350 and 2450 nm to estimate the 
contribution to the total flux emitted outside the observed range. 
The match of the normalization factors at the two ends is in all 
cases very good. We show in Fig.[5]an example of this procedure. 

The main source of uncertainty in L* comes from the un- 
certainty in the spectroscopic flux. For the objects observed in 
excellent weather conditions this matches the fotometric fluxes 
to better than a factor <1.5. It is then reasonable to assume such 
an uncertainty for all the objects in our sample, after normaliz- 
ing the spectroscopic flux to the photometry (see Sect. 12.3) . This 
would lead to an uncertainty of less than 0.2 dex in logL*. 

To convert the bolometric fluxes obtained in this way in L» 
we adopt these distances: we assu me that YSO s in th e <xOri 
region have a distance of 360 pc dBrown et al.L 1994)). those 
in TW Hya the di stances listed in | \\feinberger et al.l ( 2012b . in 
iTorres et al l (12008). and in Mamaiek (2005), and those in Lupus 
III of 200 pc (ICom erdn. 2008), as reported in Table Q] The de- 
rived stellar luminosities are listed in Table [3] These values are 
comparable with those obtained with photometric data in the lit- 



TWA14 + BT-Settl ( T =3780) 




-13 - 



Wavelength (nm) 



10000 



Fig. 5. Example of the combination of a flux-calibrated and tel- 
luric removed X-Shooter spectrum (black) with model spectra 
with the same T e ff (green), which is normalized at the red and 
blue edges of the X-Shooter spectrum and shown only outside 
the X-Shooter range. The telluric bands in the NIR are replaced 
with a linear interpolation (red). In this figure we show the ex- 
ample of TWA 14. 
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Table 4. Spectral types obtained using the method based on the spectral indices described in Sect. [3]and in Appendix IE1 



Name 


SpT" 


VlS-ind" 


H 2 O.K2 


H 2 


sH 2 O j 


sH 2 O k 


sH 2 O h1 


1, 






HP 


TWA9A 


K5 






















S0879 


K7 






















TWA6 


K7 






















TWA25 


MO 






















TWA 14 


M0.5 






















TWA13B 


Ml 


M3.3 








M1.0 


M1.2 




MO 


5 




TWA13A 


Ml 


M3.3 








M1.0 


M0.7 




MO 


7 




TWA2A 


M2 


M3.3 


M2.0 




M2.6 


M1.5 


M0.1 


M1.4 


MO 


5 




Szl22 


M2 


M3.4 


M2.2 




M2.6 


M2.7 


M2.4 


M2.5 


M2 







TWA9B 


M3 


M3.7 


M4.0 




M5.1 


M3.1 


M3.4 


M4.4 


M3 


4 




TWA15B 


M3 


M3.8 


M4.5 




M6.7 


M3.8 


M3.3 


M5.0 


M2 


1 




TWA7 


M3 


M3.8 


M2.6 




M3.5 


M4.3 


M3.3 


M3.3 


M3 


7 




TWA15A 


M3.5 


M3.8 


M3.8 




M7.2 


M4.1 


M4.0 


M5.3 


M2 


5 




Szl21 


M4 


M4.3 


M3.0 




M3.6 


M1.2 


M5.8 


M2.9 


M2 


1 




Sz94 


M4 


M3.7 


M4.3 




M4.9 


M3.7 


M3.0 


M4.9 


M4 







S0797 


M4.5 


M4.7 


M2.3 




M2.1 


M6.0 


M5.0 


M3.6 


M5 


9 




S0641 


M5 


M5.2 


M5.3 


M5.8 


M4.7 


M4.2 


M5.1 


M5.9 


M6 


6 




Par-Lup3-2 


M5 


M5.0 


M5.7 


M5.7 


M5.8 


M4.0 


M4.5 


M6.1 


M5 







S0925 


M5.5 


M5.5 


M7.3 


M6.3 


M3.9 


Ml. 5 


M5.0 


M6.2 


M7 


3 




S0999 


M5.5 


M5.4 


M3.6 


M5.8 


M3.9 


M5.7 


M5.4 


M6.1 


M6 


5 




Szl07 


M5.5 


M5.5 


M5.7 


M5.9 


M6.1 


M6.1 


M5.6 


M6.4 


M4 


5 




Par-Lup3- 1 


M6.5 


M6.4 


M6.6 


M7.1 


L0.8 


M5.7 


M8.6 


LO.O 


M5 


5 




TWA26 


M9 




M8.4 


M8.3 


LO.O 


M7.7 


M9.9 


L0.4 


M6 


8 


M9.1 


TWA29 


M9.5 




M9.4 


M8.7 


M9.3 


M8.2 


L0.8 


L0.2 


M7 


2 


M9.5 



Notes. <a) SpT derived in this work as explained in Sect. 13.11 (W Results obtained using the spectral indices in the VIS part of the spectrum, as 
explained in Sect. 13.21 All the other columns refer to the results obtained using NIR spectral indices, as explained in Appendix|B] SpT are reported 
only in the range of validity of each index. 

Table 5. Spectral indices from lRiddick et all (120071 et reference therein) adopted in our analysis for spectral type classification. 



Index 


Range of validity 


Numerator [nm] 


Denominator [nm] 


VO 7445 


M5-M8 


0.5625 (735.0-740.0) + 0.4375 (751.0-756.0) 


742.0-747.0 


V0 2 


M3-M8 


792.0-796.0 


813.0-815.0 


c81 


M2.5-M8 


811.5-816.5 


(786.5-791.5)+(849.0-854.0) 


Rl 


M2.5-M8 


802.5-813.0 


801.5-802.5 


R2 


M3-M8 


814.5-846.0 


846.0-847.0 


R3 


M2.5-M8 


(802.5-813.0)+(841.5-846.0) 


(801.5-802.5)+(846.0-847.0) 


TiO 8465 


M3-M8 


840.5-842.5 


845.5-847.5 


PC3 


M3-M8 


823.5-826.5 


754.0-758.0 



erature, with a typical difference < 0.2 dex. Hence, they are con- 
sistent with our determinations, within the errors. 

4.2. Stellar mass and age 

In Fig. [6] we show the Hertzsprung-Russell diagram (HRD) of 
our PMS stars, built using T e ff as reported in Table [3]and L„ de- 
rived in Sec 14.11 We assign M, and a ge to our PMS sta rs by 
interpolating evolutionary tracks from Baraff e et al.l (Il998l) in 
the HRD. The resulting M* are reported in Table [3] For four 
objects (Szl07, Szl21, TWA26, and TWA29) the position in 
the HRD implies an age <1 Myr, where theoretical models are 
know n to be very uncerta in, and in fact are typically not tabu- 
lated ( Baraffe et al.li l998). We estimate the mass of these objects 
extrapolating from the closest tabulated points, but we warn the 
reader that the values are affected by high uncertainty. 

Our YSOs are distributed along different isochrones; Lupus 
YSOs appear to be younger than the others (age<2 Myr), while 
crOri YSOs are distributed in isochronal ages in the range 
2<age<10 Myr; finally, TW Hya targets appear generally close 
to the 10 Myr isochrone. This is in general agreement with 
what found in the literature; indeed, Lupus has an estimated 



age of ~ 1 - 1 .5 Myr ( Hughes eTall 1 1994 IComeronl 120081) . while 
the crOri region is usually considered to be slightly older (~3 
Myr in average), and ranging from <1 Myr to several Myr 
(Zapater o Osorio etaU 120021; lOliveira et al l|2^.ForflieTW 
Hya association, the age estimates are >10 Myr (Mamaiek, 
120051: iBarrado Y Navasculi l2006t FWeTnberger et al.lboi2b . 

5. Line classification 

The spectra of our objects are characterized by the presence of 
photospheric absorption lines that depend on the SpT and in 
some cases on the age. In order to assess the PMS status of the 
objects in our sample, we check that the lithium absorption fea- 
ture at A 670.8 nm, w hich is related to the age of the YSOs (e.g. 
iMentuch etaU 120081) . is detected in all but one (Sz94) of the 
objects. We discuss in more detail the implications of the non- 
detection in Sz94 in Appendix[A] and we explain why this object 
could be anyway considered in our analysis as YSO. The values 
of the lithium equivalent width (EWui) for the other objects in 
the sample are ~ 0.5 A. A detailed analysis of this line and the 
other photospheric ab sorption lines of the objects in our sample 
will be carried out bv lStelzer et al.l (12013. in prep.l) . 
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Fig. 6. Hertzsprung-Russell Diagram of the Class III YSOs of 
this work. Data p oints are compared w ith the evolutionary tracks 
(dotted lines) by Baraf fe et al.l ([1998). Isochrones (solid lines) 
correspond to 2, 10, 30 and 100 Myr. 



In addition to these, we detect many emission lines, typically 
H, He and Ca lines, that originate in the chromosphere of these 
stars. In this work we concentrate on the emission lines charac- 
terization, as we are interested in the chromospheric activity. 

5. 1. Emission lines identification 

In order to understand the contribution of the chromospheric 
emission to the estimate of L acc through the luminosity of 
accretion-related emission lines, we first identify in our spec- 
tra the lines typically related to accretion processes in Class II 
YSOs. Here, we describe which lines we detect and report their 
fluxes and equivalent widths in Table |6]and|7] 

The most common line detected in YSOs is the Ha line 
at 656.28 nm, that is present in the spectra of all our ob- 
jects. Emission in this line has been used as a proxy for YSO 
identification and has been related to a ccretion processes (e.g. 
iMuzerolle etail 119981: iNatta et all 120041) . This line is also gen- 
erated in chromospherically active YSOs (e.g. White & Basri, 
12.0031) . Similarly, the other hydrogen recombination lines of the 
Balmer series are easily detected in almost all of our Class III 
objects up to the H12 line {A 374.9 nm). It is not easy, never- 
theless, to determine the continuum around Balmer lines beyond 
the H9 line (A 383.5 nm), and the He line (A 397 nm) is blended 
with the Call-K line. An example of a portion of the spectrum 
from H/3 to H12 is shown in Fig. [7] 




0.0 



380 



400 420 440 460 
Wavelength [nm] 



480 



Fig. 7. Portion of the spectrum, showing emission in all Balmer 
lines from H/3 up to HI 2, as well as the Call H and K lines, of the 
YSO TWA13B. The spectrum has been smoothed to a resolution 
R = 3750 at 375 nm. 



The hydrogen recombination lines of the Paschen and 
Brackett series, in particular the Pa/3 (A 1281.8 nm) and Bry 
(A 216 6 nm) lines, have been shown to be related to accre- 
tion by IMuzerolle et al.l (Il998l) . These lines have subsequently 
been used to survey star forming regions with high extinction 
(Natt a et all |2~004, 2006) in order to obtain accretion rate esti- 
mates for very low mass objects. We do not detect any of these 
lines in our Class III spectra, confirming that chromospheric ac- 
tivity is not normally detectable with these lines. 

The Calcium II emission lines at AA 393.4, 396.9 nm 
(Ca HK) and at AA 849.8, 854.2, 866.2 nm (Ca IRT) 
are related to accretion processes (e.g. iMohantv et all [2005; 
IHerczeg & Hillenbrand l2008t iRigliaco et alll2012l) , but also to 
chromospheric activity (e.g. iMontesl Q998). The Call H and K 
lines are detected in 90% of our objects. The Call IRT lines are 
detected in 1 1 out of 13 objects with SpT earlier than M4. These 
emission lines appear as a reversal in the core of the photospheric 
absorption lines. For all 1 1 objects with SpT M4 or later, the Call 
IRT lines are not detected. 

The Hel line at A 587.6 nm is also known to be 
associated with accret i on pr ocesses (IMuzerolle et al.l |l998; 
IHerczeg & Hillenbrand [2008), but in Class III YSOs is known 
to be of chromospheric origin (e.g. lEdwards et al.l 120061) . The 
line is indeed detected in 22 (92%) objects. Other Hel lines at 
AA 667.8, 706.5 , and 1083 nm are usua l ly associated with accre- 
tion processes (IMuzerolle et all 1 19981 IHerczeg & Hillenbrand 
|2008[ lEdwards et all 120061) : we detect only in Szl22 the Hel 
lines at AA 667.8, 706.5 nm, while we detect in 8 (33%) of the 
objects the Hel line at A 1083 nm. 

Finally, there is no trace of forbidden emission lines in any 
of our X-Shooter spectra, consistent with the expected absence 
of circumstellar material in Class III YSOs. 



5.2. Ha equivalent width and 10% width 

A commonly used estim ator for the ac t ivity in PMS stars is the 
EW of the Ha line (e.g. White & Basri, 2003). This is useful es- 
pecially when dealing with spectra that are not flux-calibrated 
or with narrow-band photometric data. The absolute values of 
this quantity as a function of the SpT of the objects are plotted 
in Fig. |8] and the values are reported in Table [6] We observe a 
well-known dependence of EWn a with SpT that is due to de- 
creasing continuum flux for cooler atmospheres. With respect to 
the threshold to distinguish between accreting and non-accreting 
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s 




1000 



10% Ha Width [km/s] 



Fig. 8. Ha equivalent width as a function of spectral type. The 
dashed lines represen t the boundary betwee n accretors and non 
accretors proposed by IWhite & Basril (12003b for different SpT. 



YSOs proposed by White & Basri (2003), all our targets satisfy 
the criteria of IWhite & Basril (12003) for being non-accretors . 

Another diagnostic to distinguish between accreting and 
non-accreting YSOs is the fu l l widt h of the Ha line at 10% of 
the line peak (I White & Basril 120031) . This diagnostic has been 
shown to be corre lated with M acc , but with a large dispersion 
dNatta et aUl2004l) . Fig.|9]shows the EW Ha absolute values ver- 
sus the 10% Ha width. We see that for most of our objects the 
10% Ha width is in the ~ 100-270 km/s range, and for only 
three obj ects this value i s sign ificantly above the threshold sug- 
gested by|Wti!&!asri|(|2003|) of 270 km/s (TWA6, Szl22, and 
Szl21). In Appendix lAl we discuss about these objects, and we 
explain why TWA6 and Szl21 can be considered in our analysis, 
while Szl22 should be excluded because it is probably an unre- 
solved binary. This is probably due to high values of vsinz for 
these objects, which broaden the line profile or to unresolved bi- 
narity. For BDs, the threshold to distinguish accretors from non- 
accre tors is at 10% Ha width of 200 km/s ( Javaward hana et all 
l200l . This is satisfied for all the BDs in the sample. The values 
of 10% Ha width are listed in Table 

5.3. Line luminosity 

We measure the flux of each line by estimating the continuum in 
the proximity of the line with the IDL astrolib outlier-resistant 
mean task resistanLmean. We then subtract the continuum 
from the observed flux and calculate the integral, checking that 
the whole line, including the wings, is included in the computa- 
tion. To compute Lii ne we adopt the distances reported in Table|3] 
For the Call IRT lines, where the emission appears in 
the core of the absorption feature, we subtract the contin- 
uum from our spectrum following the prescription given by 
ISoderblom et alJ (Il993l) . Using a BT-Settl synthetic spectrum 
dAllard et all 201 ll) of the same T e g smoothed at the same res- 
olution of our observed spectrum, we obtain an estimate of the 
line absorption feature that is then subtracted in order to iso- 
late the emission core of the line; finally, we integrate over the 
continuum subtracted spectrum. We report in Tables |6]and|7]the 
values obtained for the fluxes and the line EWs. 



Fig. 9. Ha equivalent width as a function of the 10% Ha width . 
The vertical dashed line rerpesents the IWhite & Basril ({2003) 
criterion for the boundary between accretors and non-accretors. 
The objects with 10% Ha width bigger than 270 km/s are, from 
right to left: Szl22, Szl21, TWA6 and TWA 13 A. 

We include in Table [6] the values of the observed Balmer 
jump, defined as the ratio between the flux at ~360 nm 
and that at ~400 nm. Typical values found in the lit- 
erature for ClassIII YSOs range between ~0. 3 and 0.5 
(IHerczeg & Hillenbrand! l2008tlRigliaco et alll2012l) . In ClassII 
YSOs, instead, t he observed Balmer j ump values are usually 
larger, up to ~6 (|Hartigan et all 1199 It: IHerczeg & Hi llenbrand, 
2008; Rig liaco et all 120121) . For the objects in our sample, this 
quantity ranges between ~0.35 and ~ 0.55 (see Table [5J, with 
the exception of Szl22. The values of the Balmer jump ratio for 
the three BDs in the sample are not reported, because the SNR 
of the UVB spectrum of these sources is too low to estimate this 
quantity. 

6. Implications for mass accretion rates 
determination 

The physical parameter that is used to estimate the accre- 
tion activity in Class II YSOs is M acc , which is derived from 
L acc and the stellar p arameters using the following relation 
(Hartman netaQll998l) : 



Rm 



GM, 0.8GM, 



(1) 



where the factor 0.8 is due to the assumption that the accretion 
flows arise from a magnetospheric radius R m ~ 5R* (IShu et all 
[1991 . 

When estimating L acc in Class II YSOs wit h the direct 
method of UV-excess fitting (e.g. IValenti et all 1 19931) . the 
contribution to the continuum excess emission due to chro- 
mospheric activity is probably negligible and it is normally 
taken into account using as a temp late for the analysis a 
Class III YSO of th e same SpT (e.g. IHerc zeg & Hillenbrand!, 
12008! iRigliaco et all l201l! 120121) . However, L acc is often de- 
rived using spectral lines luminosity and appropriate empiri- 
cal relations (see e.g. iMuzerolle et all 119981; iNatta et all l2004t 
IHerczeg & Hillenbrand! I2008L and references therein). 
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Table 6. Fluxes and equivalent widths of Balmer lines. 
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n nnan _i_ n nm i 
U.UUjU ± U.UUIj 


^n nm o 
<U.UU1Z 


^n nm i 
<U.UU1 1 


^n nnno7^ 

<u.uuuy / j 


^n nnno77 
<U.UUUo / / 


1 /IO o 

14Z.Z 


n ao 
U.oo 




— O. JJ 


A O 1 

— o.ol 


^ ^n 

— j.jy 


a on 

— 4.zy 
















n on j- n no 
U.ZU ± U.UZ 


n m i j. n nn 1 } 
U.UjI ± U.UUj 


n n 1 A! j. n nno 
U.Ulo ± U.UUZ 


n nno 1 _i_ n nnon 
u.uuy i ± u.uuzu 


n nn^n _i_ n nm 1 
U.UUjU ± U.UU1 1 


n nmA _i_ n nnnA 
U.UUjo ± U.UUUo 


n nno/i _i_ n nnn^ 
U.UUZ4 ± U.UUUj 


n nno 1 _i_ n nnn7 
U.UUZ1 ± U.UUU/ 


1/17 1 
14 / . 1 


n /i< 
U.4j 




— o. /U 


— / .yy 


7 on 
— / .oU 


J.J 1 


A TX 
— 4.Zj 


A 1 

— 4. lo 


— Z.4j 


1 

— 1. JJ 






i>ZiU / 


i q _i_ n o 
l.o ± U.z 


n A(\ _i_ n n^ 
U.4U ± U.Uj 


n o 1 _i_ n no 
U.zl ± U.Uz 


n 1 1 _i_ n no 
U.Ij ± U.Uz 


n nA/i _i_ n m n 
U.Uo4 ± U.U1U 


n n/io _i_ n nn7 
U.U4y ± U.UU/ 


n n^A _i_ n nnA 
U.Ujo ± U.UUo 


n no^ _i_ n nno 
U.UZj ± U.UUo 


O^/l Q 

Zj4.y 


n qo 
U.jo 




-11.72 


-11.77 


-9.88 


-6.91 


-4.28 


-4.23 


-2.57 


-2.07 






Par-Lup3-1 


0.22 ± 0.02 


0.024 ± 0.007 


0.016 ±0.007 


0.0073 ± 0.0060 


<0.0082 


<0.0066 


<0.0072 


<0.0055 


140.0 


1.43 




-15.78 


-16.82 


-29.99 


-9.01 














TWA26 


0.052 ± 0.008 


0.012 ± 0.002 


0.0085 ±0.0016 


0.0068 ± 0.0017 


0.0038 ± 0.0007 


<0.0020 


<0.0022 


<0.0018 


138.4 


1.55 




-9.91 


-27.07 


-66.62 


-35.39 


-97.42 












TWA29 


0.019 ± 0.003 


<0.0016 


<0.0014 


<0.0018 


<0.0021 


<0.0017 


<0.0028 


<0.0028 


114.9 


1.56 




-14.24 





















Notes. Fluxes (10 14 erg s 1 cm 2 ) are reported in the first rows for each objects with their errors, while equivalent widhts (A) in the second rows. The Ha 10% width is in km/s. The last column 
refers to the Balmer jump ratio, intened as the ratio between the flux at -360 nm to the flux at -400 nm. Upper limits are reported with <. 



Table 7. Fluxes and equivalent widths of Helium and Calcium lines. 



Name 


Hel^587.6 


1161,11083 


CaII^3t)34 




^-^11^854.2 




1 WAyA 


5.71e-15 + 3.31e- 


15 


, 1 1 f\„ 1 A 

<l.lUe-14 


9.59e-14 ± 7.44e- 


16 


4.jJe-lJ ± o.zle-lj 


OOOd l^ _i_/1 ,-1 7c |C 

z.z/e-lo ± 4.4/e-lD 


^n a 1 "3 _i_ A 1 A ^ 1 c 




-0.08 






-13.46 




n 1 a 


n ^0 
U.Jo 


n At. 
U.4o 




1.08e-15 ±6.52e- 


16 


O.oje-lo 


9.82e-15 ± 1.57e- 


16 


a n, 1 1 a _i_ 1 n a ia 
4.U4e-14 ± o.lje-lo 


1 7C a 1 A _i_ "2 O/ld 1 £ 

1. /oe-14 ± j. /4e-lo 


O OAa 1 A _1_ Q A^a 1 A 

Z.Zoe-14 ± j.oDe-lD 




-0.15 






-18.79 




n (\a 
-U.U4 


n oq 
U.Zo 


n 11 


1 WAo 


1.87e-14 ± 1.47e- 


14 


O QQa 1 A _1_ Q QOa 1 < 

z.yye-14 ± y.y/e-ij 


9.66e-14 ± 4.87e- 


15 


<j.54e-l j 


<4.oye-l j 


.0 or a 1 c 

<Z.oDe-l 5 




-0.25 




n ^n 
— U.jy 


-12.29 










1 WAZJ 


2.39e-14 ± 1.14e- 


14 


1 O/la 1 A _l_ ^ O^a 1 A 

3.oUe-14 ± o.oje-14 


2.46e-13 ± 1.60e- 


15 


^^o. 1 1 _l_ £L A 1 a. 1^ 

o.jJe-lJ + o.4le-lj 


1 Q£.n. 1 1 _i_ n /1 Oa 1 c 






-0.19 




n 1 a 
— U. ID 


-20.04 




n n 1 

U.Ul 


n 0^ 
U.Zj 


U. j J 


T\X7A 1 A 
1 WA14 


6.27e-15 ± 5.04e- 


15 


l.ole-14 ± D.Dle-lj 


3.92e-14 ± 1.13e- 


15 


t ?/1a 1 1 l 1 OOa 1 ^ 

i.j4e-io ± i.zye-ij 


1 1 /I _i_ O A a 1< 

/.jje-14 ± Z.4oe-15 


a /iOa 1 /i _i_ 1 n^d 1 ^ 
o.4Ze-14 ± l.U^e-lj 




-0.31 




n /i/i 

— U.44 


-19.84 




n no 

U.UZ 


n /i 1 

U.41 


n /io 
U.4V 


1 WAIjd 


1.66e-14 ± 7.83e- 


15 


Q 1 /I _i_ ^ 1 /I 

o.jje-14 ± j.JUe-14 


1.22e-13 ± 7.24e- 


16 




Z.lJe-lJ ± o.Uoe-lj 


O 1 <a ll _l_/1 Qla 1^ 

Z.loe-1.3 ± 4.oJe-lj 




-0.24 




n ai 

— U.Dl 


-21.91 




n ni 


n 00 
u.zz 


n ^10 
U.jZ 


TWZA 1 1 A 


4.33e-14 ± 1.21e- 


14 


<l. / je-14 


1.83e-13 ± 1.18e- 


14 


o.jUe-lJ ± 3.0/e-lj 


QQa 1 Q , 1 Al a 1/1 

z.yje-ij ± i.uie-14 


^ 1 1 1 1 n Ala 1 < 

o.loe-lj ± /.4Ze-lj 




-0.61 






-15.88 




n 1 
-U. Iz 


n 1 ^ 
U.l j 


n oa 
U.Zo 


1 WAZA 


1.61e-14 ± 5.31e- 


15 


<j.y.ie-l4 


2.10e-13 ± 7.81e- 


16 


1 11q IO _l_ 1 n^a 1/1 

l.zJe-lz ± l.Uje-14 


A fl/1a 1 1 _l_ 1 QCa 1/1 


t o/i Q i^_i_nnia 1^ 

3.o4e-ij ± y.uie-ij 




-0.16 






-21.62 




n no 
U.Uo 


n 10 


n /io 
U.4Z 


C rj 1 00 

oZIZZ 


2.93e-15 ± 6.09e- 


16 


O 1 Oa 1 ^ _1_ O 1 1 < 

Z.loe-lj ± Z.loe-lj 


1.88e-14 ± 7.91e- 


16 


-1 QA„ A f. 

<4.o4e-lo 


<o. Joe-lo 


^-O Q/1a 1 A 

<z.y4e-io 




-0.44 




n nA 
— U.UO 


-19.58 










T\X7A on 
1 WAVd 


1.96e-15 + 2.73e- 


16 


<o. /4e-l j 


1.72e-14 ± 1.44e- 


16 


1 1 ^o. 1^1 _l_ 1 T'la 1 C 


^ Q/| Q 1 A 1 1 C7 Q K 

j.o4e-l4 ± l.j/e-lD 


j.o4e-i4 ± y.jDe-io 




-0.37 






-20.82 




n no 
U.UZ 


n 1 
U.lo 


n ^0 
U.jZ 


1 WAIjd 


2.63e-15 + 5.10e- 


16 


1 ariii 1 a _i_ /i ^n^ 1 ^ 
l.JUe-14 ± 4.:>Ue-lj 


1.51e-14 ± 3.17e- 


16 


£ 1 Ha 1 /I _1_ Q 1 1 A 

o.lUe-14 ± o.lje-lo 


£Qa 1 /I _1_ 1 1 < 

z.oye-14 ± looe-ij 


O 1 1 A _i_ ^ QO a |A 

z.iye-14 ± j.yoe-io 




-0.69 




n 


-27.41 




n on 
-U.ZU 


U.Uj 


n 1 q 
u.iv 


1 WA / 


1.66e-14±2.82e- 


15 


, 1 C7. 1 A 

<l.j /e-14 


1.42e-13 ± 8.57e- 


16 


Q AH a 11 1 Q C/!: a 1 C 


2 Qfla 1 ^ _1_ 1 1 Q a 1 /I 

j.vUe-io ± i.iye-14 


j.yie-u ± o.j4e-ij 




-0.35 






-25.00 




n 1 n 
-U. 1U 


n m 
u.u/ 


n 00 

U.ZZ 


1 WA1 j A 


4.40e-15 ± 6.80e- 


16 


1 11a 1 A _l_ /I fi/1 a 1 C 

l.oZe-14 ± 4.U4e-lj 


1.79e-14 ± 8.08e- 


16 


1 1 Q a 1 A _i_ O QAq 1 £ 

/.ioe-14 ± y.JUe-iD 


j.lle-14 ± l.jye-lD 


OO/la 1/1 _l_0 f\Aa. 1 A 

Z.o4e-14 ± /.U4e-lo 




-1.05 




— U. /U 


-19.29 




n ^ 1 


n no 
-U.UZ 


n 1 a 
U.lo 


oZlZl 


<1.45e-16 




•2 ^ria 1 ^ _1_ 1 Ola 1 ^ 

J.DUe-lj ± J.ZZe-lj 


1.22e-14 ± 5.84e- 


16 


- 1 OA, 1 C 

< 1 .j\JQ- 1 J 


<z.zje-i D 


-'7 nAd 1 A 

< / .Uoe-lo 








A 1 1 

— u. 1 1 


-23.33 










Sz94 


1.12e-15 ± 1.41e- 


16 


<y.l je-lo 


5.64e-15 ± 1.53e- 


16 


< I .j /e-lo 


<y.yje-io 


<j. /Ue-lo 




-0.59 






-18.92 










ou /y / 


2.92e-17 ± 7.82e- 


18 


, noo. i no 
<j.UZe+UZ 


2.48e-16 ± 2.09e- 


17 


^O O^a 1 


, 1 11. I /I 

<l.l le-lo 


<4.oye-l / 




-0.35 






-20.54 










cn^/i 1 
i>U041 


2.12e-17 ± 3.81e- 


18 


O.yoe-l / 


1.00e-16±2.47e- 


18 


- 1 on, 1 -7 

<4.ove-l / 


<o.o /e-l / 


^1 y|7a IO 

<j.4/e-l / 




-0.65 






-16.82 










O , , ,- T O 

rar— LupJ— 2 


3.00e-16 + 7.21e- 


17 


<y.yje-lo 


1.39e-15 ±4.54e- 


17 


<y.UUe-lo 


, 1 7C 1 c 

<l.ZDe-l 5 


<j.l De-lo 


-0.46 






-9.30 










S0925 


1.98e-17 ± 3.89e- 


18 


<7.64e-17 


8.55e-17 ± 7.46e- 


18 


<4.93e-17 


<7.12e-17 


<3.54e-17 




-0.79 






-18.79 












2.70e-17 ± 8.90e- 


18 


, 1 1 Q„ 1 (1 

<i.iye-io 


2.05e-16 ± 3.06e- 


17 


<y.o(Je-l / 


<l.oje-lo 


< /.sye-i / 




-0.42 






-15.56 










Szl07 


3.08e-16 ± l.ooe- 


16 


<5.19e-16 


3.34e-15 ± l./oe- 


16 


<5.03e-16 


<1.06e-15 


<3.70e-16 




-0.57 






-25.26 










Par-Lup3-1 


6.02e-17 ±5.72e- 
-2.35 


17 


<3.71e-16 


<4.40e-17 




<1.97e-16 


<2.86e-16 


<1.58e-16 


TWA26 


1.89e-17 ±9.96e- 


18 


<4.28e-16 


6.20e-17 ± 1.82e- 


17 


<1.10e-16 


<1.69e-16 


<8.50e-17 




-1.88 






-56.20 










TWA29 


<2.27e-17 




<9.49e-17 


<2.02e-17 




<2.59e-17 


<4.37e-17 


<2.73e-17 



Notes. Fluxes are reported in [erg s 1 cm 2 ] in the first rows, Equivalent Widths (EW) in [A] are reported in the second rows. Upper limits are shown with <. 



Manara et al.: Impact of chromospheric emission on accretion rate estimate of young stellar objects 



Several of the lines normally used for these studies are in- 
fluenced by chromospheric activity (Sect. [5), which should be 
estimated and subtracted from the line emission before com- 
puting L acc . This procedure is not trivial, as it is difficult to 
properly disentangle the two emission processes in each ob- 
ject and each line. All the L acc - Lii nes relations are always 
based on the non-corrected v alue s of Ly ins ,,. This correc tion, as 
shown in llnglebv et al.l d201 lb and Rig liaco et al.l d2012l) . can be 
important in objects with low L acc and in particular in VLM 
stars. Therefore, the chromospheric emission acts as a system- 
atic "noise" which affects the L acc measurements. 

In the following, we will characterize this effect using the 
Lii n e derived in Sect. [5] and th e most recent L acc - Lr me re lations 
for Class II YSOs derived bv lAlcala et"ai1(l2013. in prep.l) . These 
re lations have been deri ved using the same method as described 
in Ri gliaco et al.l d2012h . where L acc is obtained from the contin- 
uum excess emission and then compared with the L wof the sev- 
eral emission l ine diagnostics in every object. In lAlcala et al.l 
in prep.l) the sample is composed of 36 YSOs located 
in the Lupus-I and Lupus-III clouds, togeth er with 8 additiona l 
YSOs located in the cr-Ori region from Riglia ccTet al.l (1201 2b . 
These relations are not quantitatively different from those in the 
literature, but have significantly smaller uncertainties and have 
been derived for stars with similar properties than the Class III 
analyzed here. For each Class III object, we compute L acc from 
a number of different lines, as described in the following. This 
provides a measurement of the "noise" introduced in the deter- 
mination of L acc from line luminosities in Class II. It represents a 
typical threshold for the determination of L acc in Class II objects 
by chromospheric activity, assuming that this is approximately 
the same in the two different classes of objects, with and without 
ongoing accretion; we define it in the following as L acc no i Se . 

6. 1 . Accretion luminosity noise 

Alcala et afj (I2013. in prep.1) are using a sample of Class II YSOs 
in the Lupus star forming region observed with X-Shooter to 
refine the L acc - Lii ne relations. We adopt their relations to es- 
timate the L acc no i se for our Class III YSOs, using in particular 
the Ha, U/3, Hy, U6, H8, H9, H10, Hll, Hel (X1587.6, and 
1083 nm), Call (,1393 nm), Call (,1849.8 nm), Call (A854.2 nm) 
and Call (/1866.2 nm) lines. Moreover, we use the relation from 
iNatta et all (|2004|) between M acc and the 10% Ha width to esti- 



mate L; 



acc, noise 



from this indicator. 



We show in Fig. [TOlfTTI the values of L acc , no ise for every ob- 
ject obtained using the different indicators. The uncertainties on 
these values are dominated by the errors in the relations between 
Lii ne and L acc . Upper limits for undetected lines are shown with 
triangles and indicate the 3<x upper limits. Each Balmer line 
leads to values of L acc ,noise that are always in agreement with 
the other Balmer lines by less than ~0.2 dex, and similarly the 
HeLj587.6 line almost in all cases. Solid lines represent the arith- 
metic mean of the Balmer and HeL;587.6 L acc no i se values, which 
is also reported in Table [3] while the dashed lines show the one 
sigma standard deviation. The HeI,uo83 an d the Call lines, in- 
stead, are in various cases not in agreement with the result ob- 
tained using the Balmer and Hel^587 6 lines, with differences even 
larger than 0.6 dex. It should be considered that the exact value 
of the Call IRT lines luminosity is subject to many uncertainties, 
due to the complicated procedure to estimate the excess luminos- 
ity (see Sect. 1531). 

The Paschen and Brackett HI emission lines are not de- 
tected in our spectra, as pointed out in Sect. 15.11 We report 
in Fig. [TOffm the 3<x upper limits for L acc no ise obtained us- 
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Fig. 12. Mean values of log(L acc n0 i se /Lo) obtained with different 
accretion diagnostics as a function of T e ff . Error bars represent 
the standard deviation around the mean log(L acCin oise/L ). These 
data should be intended as the noise in the values of L acc due to 
chromospheric emission flux. 



ing the Pa/3 and Bry line l uminosities and the relations from 
Alc ala et al.ld20rT in prep.1) . These values are always below the 
mean L acc ,noise value obtained with the Balmer and HeF^.e 
lines. This implies that those lines are less sensitive to chromo- 
spheric activity than the Balmer lines. 

The 10% Ha width is the accretion indicator that leads to 
values of L acc 



e more discrepant from the mean (blue circles 
in Fig. ITOlfTTb . This clearly does not follow in more than 50% of 
the cases the results obtained using the other indicators. This is 
not surprising as the Ha width is mainly a kinematics measure- 
ment, contrary to Lii ne measurements. It is to be expected that 
the application of a method calibrated for accretion processes to 
chromospheric activity would result in inconsistencies. We know 
that the broadening of the Ha line and the other accretion-related 
lines is due to the high-velocity infall of material in the accretion 
flows, while the intensity of the emission lines is due to emission 
from high temperature region. The latter can be either accretion 
shocks on the stellar surface or chromospheric emission. The 
fact that in our sample of non-accreting objects relations con- 
verting Liine to L acc lead to similar results when using line fluxes, 
while the result is quite different when using line broadening 
seems to confirm that the Ha 10% width we detect is only due 
to thermal broadening in the chromosphere of these stars and 
not to gas flow kinematics associated with the accretion onto the 
central object. 

For all objects, the mean L accn oise value is always below 
~ 1O~ 3 L and this value decreases monotonically with the SpT 
In Fig. [T2]these mean values of log(L accn oise/L ) obtained with 
the Balmer and HeFi587.6 lines are plotted as a function of the 
T e ff of the objects. The error bars on the plot represent the stan- 
dard deviation of the derived values of L acc no i Sf 0- These values 
should be intended as the noise in the L acc values arising from 
the chromospheric activity. In Fig.QjJwe show the mean values 
of the logarithmic ratio L acc no i se /L* obtained using the Balmer 



For the object TWA29, where only the Ho- line is detected, we re- 
port with a dashed line the error on the estimate of L m 
line in Fig. 1 121131 
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Fig. 13. Mean values of log(L acc n0 i se /L») obtained with different 
accretion diagnostics as a function of logT e g\ The dashed line 
is the best fit to the data, whose analytical form is reported in 
Eq. (0. Two objects (Szl22 and TWA9A) are excluded from the 
fit (empty symbols), as explained in the text. 



-9.0 



-9.5 



-11.0 



-12.0, 




-0.8 



-0.6 -0.4 -0.2 

log(M/M sun ) 



0.0 



0.2 



acc.noise 



as a function of logM,, with values 
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Fig. 14. logM. 



Bara ffe et al.l £l998) and the values of L acc no i se /L» derived from 
the fit in Eq. (0 at any T e ^. Results using the 1 Myr isochrone 
are reported with filled circles, those using the 3 Myr isochrone 
with filled triangles, and those using the 10 Myr isochrone with 
filled squares. 



and Hel,i587.6 lines as a function of the T e ff. Contrary to Fig. [12] 
the quantity L a cc,noise/L, is unbiased by uncertainties on distance 
values or by different stellar ages, leading to smaller spreads. 
We see that from the K7 objects down to the BDs, the values of 
log(L acc no i se /L») decrease with the T e ff of the YSOs. Fitting with 
a power-law the L aC c,noise- T e ff relation, using only the objects in 
the range K7-M9.5 and excluding Szl22 (see Appendix [Al for 
detail), we obtain the following analytical relation (dashed line 
in Fig. [13]): 

log(L a cc,„oise/i») = (6.17 ± 0.53) ■ log r eff - (24.54 ± 1.88) (2) 

The only clear deviation from the general trend, apart 
from Szl22, is the K5 YSO TWA9A, which shows a value 
of log(L acc no i se /L*) lower by ~0.6 dex with respect to what 
should be expected by the extrapolation of the previous relation. 
Unfortunately, our sample is too small and we do not have other 
objects with earlier SpT to verify whether this low value is ac- 
tually a different trend due to different chromospheric activity 
for earlier SpT YSOs or if the source is peculiar. There are also 
signatures of different chromospheric activity intensity among 
objects with the same SpT and located in the same region; for 
example, the two TW Hya Ml YSOs, which are two compo- 
nents of a binary system, thus coeval objects, have a spread in 
logCLaccnoise/L,) of -0.5 dex. 

Values of Lr me and L acc in Class II YSOs that are of the order 
of those estimated in this work should be considered very care- 
fully, as the chromospheric activity could be an important factor 
of the excess luminosity in the line and could produce mislead- 
ing results. 

6.2. Mass accretion rate noise 

In this section, we determine what the typical "chromospheric 
noise" on M acc would be when derived from indirect methods if 
the chromospheric emission is not subtracted before computing 
Liine- We will refer to this quantity as M acc no i se . 

The procedure used is the fo llowing. We select three 
isochrones (1,3, and 10 Myr) from Baraffe et al] d!998l) mod- 
els and nine different YSOs masses (0.1 1, 0.20, 0.35, 0.40, 0.60, 



Table 8. Values of logM acc no ; se at different M* and ages. 



M, 




Age [Myr] 




[M ] 


1 


3 


10 


0.11 


-10.59 


-11.08 


-11.62 


0.20 


-10.15 


-10.80 


-11.26 


0.35 


-9.96 


-10.41 


-10.94 


0.40 


-9.90 


-10.32 


-10.83 


0.60 


-9.66 


-10.05 


-10.51 


0.75 


-9.50 


-9.90 


-10.33 


0.85 


-9.40 


-9.80 


-10.20 


1.00 


-9.28 


-9.67 


-9.97 


1.10 


-9.18 


-9.52 


-9.80 



Notes. Values of log M acc noise at different M , and different ages obtained 
using isochrones from Bar affe et al.1 JT998"> and the procedure explained 
in Sect.l6T2l 



0.75, 0.85, 1.00, and 1.10 M ), which correspond always to 
r e ff in the range 2500-4000 K, where our results are applicable. 
Then, for each r e jf we derive L acc ,noise/£» using the fit reported in 
Eq. (0. Finally, we use Eq. CD, adopting the proper R„, M„ and 
L» at any age from the Baraffe et al. (1998) tracks, in order to 
determine the typical M acc no i se at different ages as a function of 
M*. The results are shown in Fig. [14] where a strong correlation 
between the two parameters is evident, with increasing M acc no i se 
with M». At the same time, the variation of M acc , no ; se with age at 
any given M* is large, up to 0.5 dex for differences of 2 Myr at 
the H-burning limit, but decreases with increasing M,. Similar 
results are obtained when using other evolutionary models. 

We derive a limit on the detectable M acc of ~ 6.6-10~ 10 M Q /yr 
for solar-mass, young (1 Myr) objects, decreasing to 2.5 ■ 10~ 12 
M /yr for low-mass, older (10 Myr) objects. We report these 
^acc.noise values for the three isochrones analyzed in Table [8] 

6.2.1. Comparison with literature data 

As discussed, M acc no j se is a lower limit to the values of M acc that 
can be derived from the luminosity of lines emitted by stellar 
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Fig. 15. logM a[x a s a function of logM* for Class II objects l ocated in different sta r forming regions. Data for pOph a r e from 
Natt a et al.l (1200 6). corrected for new distance : 



measurem ents by Rigliaco et al] d201 lb: ONC po i nts are from lManara et al] (1201 2|) . 
data for L1630N and L1641 are from Fan g et aD (120091) . and those for Chall from Biazzo et al.1 (120121) . M acc values obtained with 
indirect methods are shown with coloured filled points, while measurements based on the direct U -band excess method are shown 
with grey empty points as a reference. Downwards triangles refer to upper limits. The thick red solid line is the lower limit to the 
measurements of M acc set by chromospheric activity in the line emission. We use the values for the correct isochrone according to 
the mean value of the age for each region, as reported on the plot. 



chromospheres. In this section, we compare this chromospheric 
limit with estimates of M acc based on Lii ne from the literature. 
We show in Fig. Q3] values of M acc as a function of M* ob- 
tained with optical or NIR emission lines as accretion diagnostic 
for objects located in four star forming regions with different 
ages: p-Ophiucus, the Orion Nebula Cluster, the L1630N and 
LI 641 regions, and the Chameleon II region. We overplot with 
red solid lines the locii of M acc ,noise obtained using the proper 
isochrone for the mean age of each region. In this way we do 
not address possible differences in M acc due to age spreads in 
the se regions, but t hese spreads, if present, are <l-2 Myr (see 
e.g. lReggiani et all 1201 lb . The effect of a similar age spread on 
the M a cc,noise threshold would be of ~0.5 dex for low-mass ob- 
jects and of ~0.3 dex for solar-mass stars (see Fig. [14}. 

We consider in this analysis M acc values obtained using op- 
tical emission lines. In Fig. [15] we sh ow the values for the -1 
Myr old regions L1630N and L1641 dFang et all 120091) . where 
the Ha, Hp 1 , and Hel lines were used, those for the ~ 2-3 Myr 
old Orion Nebula Cluster dManara et al.l 1201 2l) obtained with 
photometric narrow-band Ha luminosity estimate^ and those 
obtained with the H a, H/j, and Hel line s for the ~4 Myr old 
Chameleon II region (Bi azzo et all l2012h . In all these cases, the 
vast majority of datapoints are found, as expected, well above 
the M acc ,noise threshold, confirming that chromospheric contri- 

4 We report also the v alues of M acc obtained through {/-band excess 
by iManara et alj ( 120121) . which are shown with empty symbols as a 
comparison. 



bution to the line emission is negligible for strongly accreting 
YSOs. Nevertheless, there are a few (11) objects in L1630N and 
LI 641, where measured values of M acc are smaller than the ex- 
pected M acc , no j se . One possibility to explain this result is variable 
accretion in these objects, which has been found to vary as much 
as 0.4 dex over time-scales of 1 year (ICostigan et aUl2012l) . Still, 
this does not explain why only lower mass object happen to be 
below the threshold. In any case, these points should be con- 
sidered with caution, as the measured Lii ne could be completely 
due to chromospheric emission, leading to erroneous estimates 
of M acc . 

For the ~1 Myr old p Ophiucu s region we consid er the Lij ne 
derived through Pa/? and Bry lines (Na tta et alll2p06[). corrected 
for a more recent estimate of the distance (see iRigliaco et all 
1201 lh . As we noted in Sect. 15. 1 1 and 16. 11 Pa/? and Bry lines are 
not detected in our Class III YSOs spectra. In Fig. [15] all the de- 
tections (blue filled circles) are located well above the M acc no i S e 
locus, while the upper limits (blue open triangles) are distributed 
also at the edge of the M acc no i se threshold. This confirms the va- 
lidity of these NIR lines as good tracers of accretion and the fact 
that they are most likely less subject to chromospheric noise than 
the Balmer and Hel^s^g lines. 



7. Conclusion 

In this paper, we presented the analysis of 24 diskless, hence 
non-accreting Class III YSOs, observed with the broad-band, 
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medium-resolution, high-sensitivity VLT/X-Shooter spectro- 
graph. The targets are located in three nearby star forming re- 
gions (Lupus III, <x Ori and TW Hya) and have SpT in the range 
from K5 to M9.5. We have checked the SpT classifications, us- 
ing both spectral indices and broad molecular bands. Moreover, 
using the flux calibrated spectra, we have derived the stellar lu- 
minosity. Then, we have analyzed the emission lines related to 
accretion processes in accreting objects that are present in these 
spectra and, from their luminosities, we have studied the impli- 
cations of chromospheric activity for M acc determination in ac- 
creting (Class II) objects. This has been done by deriving the 
parameter L acc no j se , which is the systematic "noise" introduced 
by chromospheric emission in the measurements of L acc from 
line emission in Class II objects. For this analysis we assumed a 
similar chromospheric activity in the two classes of objects. 
Our main conclusions are: 

1 . All hydrogen recombination emission lines of the Balmer se- 
ries are detected in our sample of Class III YSOs spectra 
when the S/N is high enough. On the opposite, Paschen and 
Brackett series lines are not detected in emission, and they 
are significantly weaker, when compared to Balmer lines, 
then in Class II objects. The chromospheric "noise" in these 
lines is lower than in the optical lines (see Fig. fTOlfTTI i. and 
they are very good tracers of accretion in low-mass, low ac- 
cretion rates objects. 

2. Using Balmer and HeLj587.6 lines and the calibrated relations 
between Li; ne and L acc from the literature, we derive L acc ,noise 
values always in good agreement among all the lines. 

3. Calcium emission lines in the NIR spectral range (AA 849.8, 
854.2, 866.2 nm) are detected in 1 1 objects (45% of the sam- 
ple) superposed on the photospheric absorption lines. This 
results in a more complicated line flux measurement than for 
other lines. Their behavior with respect to the hydrogen line 
luminosities is different; in particular, the values of L acc ,noise 
obtained using these lines are often not in good agreement 
with those obtained using Balmer and HeL;587.6 lines. 

4. The mean values of L acc ,noise for the objects in our sample are 
smaller than ~1O~ 3 L and have a clear dependence with T e ff 
for K7-M9.5 objects. Therefore, L acc of this order or smaller 
measured in Class II objects using line luminosity as a proxy 
of accretion may be significantly overestimated if the chro- 
mospheric contribution to the line luminosity is not taken 
into account. 

5. Our results show that the "noise" due to chromospheric ac- 
tivity on the estimate of M acc in Class II YSOs obtained using 
secondary indicators for accretion has a strong dependence 
with M* and age. Typical values of log(M accno j se ) forM-type 

YSOs are in the range from 9.2 for solar-mass young (1 

Myr) objects to -11.6 M /yr for low-mass, older (10 Myr) 
objects. Therefore, derived accretion rates below this thresh- 
old should be treated with caution as the line emission may 
be dominated by chromospheric activity. 
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Appendix A: Comments on individual objects 

AA. Sz94 

Sz94 was discovered as an Ha emitting star in the survey by 
Schwarz (1977). Since then, it has been considered as a PMS 



star in reviews and inve s tigations (|Krautteri 19921: iHughes et all 



2011), but nothing was 



Il994t IComeronL |2008; Mort ier et all 
mentioned about the presence of the lithium absorption line at 
670.8 nm. Based on its spectral energy di stribution, S z 94 ha s 
been classified as a Class III IR YSO by iMerin et alj ((2008). 
For all these reasons, we included the star in our programme 
as a Class III template of M4 SpT. Notwithstanding the high 
quality of the X-Shooter data in terms of S/N and resolution, 
the Li I A 670.8 nm line is not present in the spectrum, and we 
determine an upper limit of EWu < 0.1 A. The question then 
rise whether lithium may have already been depleted in the star. 
However, lithium is significantly depleted by large factors only 
after several tens of Myr, inconsistently with the average age of 
a few Myr of the Lupus members. Analysis on the radial veloc- 
ity of this object show t hat its value is in the range of typical 
values for Lupus sources (IStelzer et all 12013. in prepl) . We con- 
sider here the source as a PMS given its position in the HRD, 
the Ha emission and the radial velocity measurement, but more 
analysis should be carried out to confirm its PMS status. 



A.2. Sz122 



Szl22 is a Class III object classified with Spitzer (IMerin et all 
2008). All the lines of this object appear very broadened. 
Nevertheless, the spectrum cannot be fitted with synthetic spec- 
tra broadened at reasonable values of vsini, meaning that this is 
not a single fast rotator. We think, therefore, that this object is a 
binary system. This could explain the very faint Lil line (EWu < 
0.25 A) and the lack of emission in the Call IRT absorption fea- 
tures, that is usually found in other early M-type objects of our 
sample. Even the presence of a broad Ho- line (10% width > 600 
km/s) can be explained by the presence of two objects. We do not 
consider this object in the analysis of the implications of chro- 
mospheric activity on M acc measurements, because the effect of 
its binarity cannot be accounted for. 



A.3. Sz121 

Similarly to Szl22, the lines of this object appear very broad- 
ened. Als o in this case, the object is classified with Spitzer as a 
Class III dMerfn et alll2008l) . We can fit the spectrum with a syn- 
thetic spectrum of the same T e ff, logg=4.0 and vsini = 70 km/s. 
This value of vsini is rather high for a typical M4-type YSO. 
If true, it would be an extreme ultra-fast rotator. The large 10% 
Ha width (~380 km/s) is then due due to its very high rotational 
velocity. Another possibility, is that this object may also be an 
unresolved spectroscopic binary. 



AA. TWA6 

This object has bee n discovered and cl assified as a young non- 
accreting object by IWebb et all |j999). It is a member of the 
TWA associ a tion an d it is a known fast rotator (vsini=12 km/s, 
ISkellv et all (|2008)). We measure a 10% Ha width of 362 
km/s. This is larger than the threshold t o distinguish a c cretin g 
and non-accreting objects proposed by White & Basril (2003). 
Nevertheless, this object can be considered a Class III YSO for 
the small EWh„, the small Balmer jump and, in particular, the 
IR classification. For these reasons we consider this object in our 
analysis. 
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Appendix B: NIR spectral indices 

Spectral indices in the NIR part of the spectrum are particularly 
useful to classify VLM stars and BDs, which emit most of their 
radiation in this spectral region. Therefore, various analyses to 
calibrate reliable NIR indices to classify late M-type objects and 
BDs have been carried out in the past (e.g. iKirkpatrick et all 
[1991 . With our sample, we are able to verify the validity of var- 
ious NIR indices for M-type YSOs, comparing the SpT obtained 
through optical spectroscopy (Sect.|3]l. We considerin this analy- 
sis different NIR indices which have been calibrated using either 
dwarfs or yo ung stars (subgiants) for ob jects with SpT M or later. 
In particular, Roias-Avala et aTl (l2012l) calibrated the H 2 - K2 
index using a sample of M dwarfs, and d erived a relation valid 
for the whole M class: lAllers et al.l(l2007l) calibrated on a sample 
of young BDs and dwarfs the gravity independent spectral index 
H 2 Q, which is va lid for objects w ith SpT in the range M5-L0; 
iTesti et aT] d200 ll) and lTestil d2009) proposed various spectral in- 
dices (sH 2 O y , sH 2 0*, sKJ, sHJ, sH 2 O m , sH 2 O m , lj, l H , and 
Ik) to classify M-, L-, and T-t ype BDs, ca librating those on a 
large sample of dwarfs; finally. IScholz et all (1201 2|) used a sam- 
ple of VLM YSOs to calibrate the HP index, which is valid for 
objects with SpT in the range M7-M9.5. We report those indices 
in Table IB~T1 

Fig. lB.ll shows the results of this analysis: the SpTs obtained 
in Sect. [3] are reported as a function of each spectral index for 
our objects (red crosses when SpTs are in the range of validity 
of an ind ex and red circles whe n not). For the spectral indices 
derived in Tes ti et all d2001l) and lTestil d2009l) . we report also the 
data for the sample of L- and M-type dwarfs they adopted in the 
analysis (green crosses). We stress that the latter sample should 
be considered with cautions, as gravity dependent spectral in- 
dices calibrated using samples of dwarfs may not be reliable for 
YSOs, which are subgiants. The published relations between the 
SpT and the spectral indices in the range of validity of each index 
are shown as black dashed linesS 

We see that the indices sH 2 O m , sKJ, sHJ, and 1^ cannot be 
used to classify YSOs of M-type class. On the opposite, we find 
a good correlation with the indices sH 2 0^, sH 2 O y , sH 2 O m , I/, 
and Ifj. For these indices we decide to fit altogether our objects 
and those from the literature in all cases, using either a linear or 
a second degree polynomial fit. We show in Fig. IB. li the best fit 
obtained, with a black solid line in case of linear fit and a red 
solid line when the best fit is polynomial. Using these relations, 
we derive the SpT for each of our YSOs, and we report these 
results in Table |4] We thus propose new SpT-spectral index re- 
lations for th e following five indices from lTesti et al. I d200l and 
iTestil (120091) . The first three are valid for objects with SpT later 
and equal to Ml: 

SpT - code = 0.55 + 2.48 • sH 2 K (B.l) 

SpT - code = 0.38 + 2.89 • sH 2 O m - 1.56 • (sH 2 O m ) 2 (B.2) 

SpT- code = 1.11 + 3.92 ■ I H - 5.35 ■ (I H ) 2 (B.3) 

The folloewing two indices are valid for objects with SpT 
later and equal to M2: 

SpT - code = 0.70 + 2.90 • sH 2 J - 1 .78 • (sH 2 J f (B.4) 

SpT - code = 1.83 + 1.08 • Ij - 2.59 • (Ij) 2 (B.5) 

where the SpT are coded in the following way: M0 = 0.0, M9 
0.9. L5 = 1 .5 and a variation of 0. 1 corresponds to a step of one 

5 No published relations are available for the indices I j and l K 



subclass. For these spectral indices we conclude that results for 
objects with SpT in the nominal range of validity of each index 
are reliable within a typical uncertainty of ~ one subclass. 

A good correlation is also found for the spectral index 
H 2 0-K2; using the analytical relation between SpT and spectral 
index from the literature, we obtain for 13 out of 22 objects with 
M SpT results that are compatible within one sub-class with the 
correct SpT (see Table|4]i. The differences can be due to the non 
perfect telluric removal in the first interval of interest of this in- 
dex. Moreover, we confirm that the H 2 index is valid for YSOs 
with SpT in the range M5-M9.5, finding an agreement within 
one sub-class for all our objects (see Table|4]i. Regarding the HP 
index, we confirm that it is not valid for YSOs with SpT earlier 
than M7, and we observe that the SpT obtained with this index 
confirm those from literature for our two later SpT objects (see 
Tableg}. 
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Spectral index values 

Fig. B.l. Spectral type of the objects as a function of the spectral index values obta ined with differen t NIR indices (Table iBTTT i. Red 
symbols are values from this work, while green crosses are from lTesti et al.l d2001) and lTestil d2009t) . Red stars are used for objects 
with SpT in the range of validity of the index, while red pentagons for those not in the range of validity of the index. Dashed lines 
are the best fit from the reference cited in Table IB. fl for each index, when available. Black and red solid lines are the best fits (black 
= linear, red = polynomial) from this work, considering both our sample and the values from the literature. 



Table B.l. NIR spectral indices analysed in Appendix iBl 



Index 


Range of validity 


Numerator [nm] 


Denominator [nm] 


Reference 


H 2 0-K2 


M0-M9 


(2070-2090) / (2235-2255) 


(2235-2255) / (2360-2380) 


Roias-Avala et al. (2012) 


h 


M0-T9 


(1090-1130) + (1330-1350) 


2-(1265-1305) 


Testi (2009) 


Ih 


M0-T9 


(1440-1480) + (1760-1800) 


2-(1560-1600) 


Testi (2009) 


Ik 


M0-T9 


(1960-1990) + (2350-2390) 


2-(2120-2160) 


Testi (2009) 


H 2 


M5-L5 


1550-1560 


1492-1502 


Allers et al. (2007) 


sHJ 


L0-L9 


(1265-1305) - (1600-1700) 


0.5-[(1265-1305) + (1600-1700)] 


Testi etal. (2001) 


sKI 


L0-L9 


(1265-1305) -(2120-2160) 


0.5-[(1265-1305) + (2120-2160)] 


Testi etal. (2001) 


sH 2 cy 


L0-L9 


(1265-1305) - (1090-1130) 


0.5-[(1265-1305) + (1090-1130) ] 


Testi etal. (2001) 


sH,O h1 


L0-L9 


(1600-1700) - (1450-1480) 


0.5-[(1600-1700) + (1450-1480)] 


Testi etal. (2001) 


sH 2 0" 2 


L0-L9 


(1600-1700) -(1770-1810) 


0.5-[(1600-1700) + (1770-1810)] 


Testi etal. (2001) 


sITO* 


L0-L9 


(2120-2160) - (1960-1990) 


0.5-[(2120-2160) + (1960-1990)] 


Testi et al. (2001) 


HIP 


M7-M9.5 


(1675-1685) 


(1495-1505) 


Scholz etal. (2012) 
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Appendix C: On-line material 

C.1. NIR spectra 
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Fig. C.l. Spectra of Class III YSOs with spectral type earlier than M3 in the NIR arm. All the spectra are normalized at 1700 nm 
and offset in the vertical direction by 0.5 for clarity. The spectra are also smoothed to the resolution of 2000 at 2000 nm to make 
easier the identification of the features. 
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Fig. C.2. Same as Fig. IC.21 but for spectral types between M3 and M5. 
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Fig. C.3. Same as Fig. IC.3L but for spectral types later than M5. 
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C.2. UVB spectra 
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Fig. C.4. Spectra of Class III YSOs with spectral type earlier than M2 in the UVB arm. All the spectra are normalized at 450 nm 
and offset in the vertical direction for clarity. The spectra are also smoothed to the resolution of 1500 at 400 nm to make easier the 
identification of the features. 
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Fig. C.5. Same as Fig. IC.5L but for spectral types between M2 and M4. 
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